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We calculate the nuclear enhancement in the angular dis- 
tribution of Drell-Yan pairs produced in proton-nucleus reac- 
tions. Nuclear effects are encoded in universal twist-4 parton 
correlation functions. We find that the Lam- Tung relation for 
the angular coefficients of the lepton-pair distribution holds 
for the double-hard, but not for the soft-hard contribution. 
We also predict that nuclear enhancement effects at RH1C 
energies can be large. 

The factorization theorems of QCD Q permit a rigor- 
ous treatment of hard processes by absorbing soft con- 
tributions to scattering matrix elements into universal, 
process-independent distribution functions. The mea- 
surement of these distribution functions in one process 
facilitates the prediction of cross sections for other pro- 
cesses. An example is the Drell-Yan process where two 
hadrons scatter to produce a lepton pair of high invari- 
ant mass Q 2 . If Q 2 is large enough, the factorization 
theorems state that the inclusive Drell-Yan (DY) cross 
section can be expressed as the product of a calculable 
hard part and universal twist-2 parton distribution func- 
tions that can be measured, for instance, in deep- inelastic 
lepton-nucleon scattering. 

However, the structure of strong interaction dynam- 
ics is much richer than can be described solely in terms 
of twist-2 structure functions, which are two-point cor- 
relators and have a probabilistic interpretation in terms 
of partons. In reactions involving only isolated hadrons, 
hard processes are usually dominated by a single hard 
scattering, and correlation functions involving more than 
two fields are normally suppressed by powers of the large 
scale. This is different when reactions are considered 
where individual hadrons interact with nuclei, or one nu- 
cleus interacts with another. In that projectile 
parton can scatter twice within the same nucleus of mass 
number A, and the double-scattering cross section is en- 
hanced by a factor A 1 / 3 , compensating for some of the 
inherent suppression of multi-parton interactions |^,^| . 

Although we will focus here on proton-nucleus inter- 
actions, our work has been primarily motivated by the 
fact that higher- twist effects should be most prominent in 
nucleus-nucleus collisions at high energy, as they will be 
routinely studied at the Relativistic Heavy Ion Collider 
(RHIC). Many experiments at RHIC will focus on a Q 2 
range where higher-twist contributions could be signifi- 
cant g . Higher-twist effects will be especially important 
in the production of minijets, which are thought to pro- 
vide the main mechanism for the initial energy deposition 



in the central rapidity region J5|,^| . 

The determination of multi-parton correlation func- 
tions in nuclei is also an important step towards the ap- 
plication of QCD to nuclear physics, in general. The en- 
hancement of higher-twist effects in reactions involving 
nuclei facilitates the extraction of these correlation func- 
tions, which encode essential aspects of the difference in 
the quark-gluon structure between isolated hadrons and 
nuclei. However nuclear enhancement could point to a 
general problem, as it might signal an early breakdown 
of perturbative QCD as applied to high-energy collisions 
involving nuclei. 

In the following we will discuss the Drell-Yan process 
initiated by a single hadron h scattering off a nuclear 
target with mass number A: 

h{P 2 ) + A{P 1 A)^l+l-+X (1) 

At leading-twist level, a parton from the beam hadron 
reacts with a single parton from the nucleus. Figure 0(a) 
shows a typical twist-2 diagram describing production of 
the lepton pair with a high transverse momentum relative 
to the beam axis. It was already pointed out some time 
ago by Guo JjJ that higher-twist effects in this process are 
enhanced in this kinematical region (see also || ) . As the 
parton of the beam hadron travels through the nucleus it 
can scatter off partons from other nucleons, generating a 
strong dependence on the size of the nucleus. Such a de- 
pendence has been observed in several experiments 
at a scale hard enough to resolve the partonic substruc- 
ture and thus to allow for a treatment in the framework 
of perturbative QCD. In Fig. |l|(b) we give an example 
for double scattering. 

To extract higher-twist correlation functions from ex- 
perimental data it is worthwhile to consider observables 
for which such contributions are large. Ideally, one is 
looking for an observable with a vanishing contribution 
from single scattering. Such quantities can be found in 
the angular distribution of DY pairs. Therefore, we here 
extend the calculation of Ref. |Q to the differential DY 
cross section da/dQ 2 dq\ dydfi, keeping the full angular 
dependence. In addition, we apply our results to p + A 
collisions at RHIC. Some time ago, Lam and Tung de- 
rived a relation similar to the Callan-Gross relation 
in deep-inelastic scattering, which states that the longi- 
tudinal helicity amplitude Wl for the virtual photon in 
the DY process is exactly twice as large as the double 
spin- flip amplitude Waa- This is a leading- twist predic- 
tion valid up to order a s . 
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FIG. 1. (a) Schematic single- and (b) double-scattering di- 
agrams for a hadron h colliding with a nucleus A and produc- 
ing a lepton pair. One additional unobserved parton 
is radiated with high transverse momentum q±. 

Due to the complex structure of multi-parton correla- 
tions a larger number of these correlation functions exist 
than at leading twist. It is useful to consider observables 
which receive contributions from only a limited number 
of correlation functions. We find that the number of con- 
tributing correlation functions can be minimized by using 
different lepton-pair center-of-mass frames (see Ref. [[l2| ) . 
For the DY process, the so-called Gottfried- Jackson (GJ) 
frame turns out to be the most convenient choice. How- 
ever, comparing the angular correlations in the GJ and 
the Collins-Soper (CS) frames also gives valuable insight 
into the nature of double-scattering contributions Jl9| . 

The differential DY cross section is a classical example 
for a two-scale process in QCD, because both Q 2 and g 2 , 
the transverse momentum of the photon, are detected. 
In order to avoid the occurrence of large logarithms of 
the type log 2 (Q 2 /<7 2 ), which are common in this case, 
we require that both scales are of the same order. This 
condition also allows us to neglect interference terms that 
would introduce additional twist-4 matrix elements g] . 

The differential cross section for DY production is 
given by 
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where S is the center-of-mass energy squared, L^ v and 
W^ v denote the leptonic and hadronic tensors, and the 
angles 9 and <$> in df2 refer to the polar and the azimuthal 
decay angles of the lepton pair in the rest frame of the 
virtual photon. The remaining freedom to choose these 
angles gives rise to the different choices of, e.g., the GJ 
and CS frames. The hadronic tensor is given by 



WV = / d A xe^ x {P x P 2 \j lx {x)j v {Q) \P X P 2 ) 



(3) 



where P\ and P 2 are the four-momenta of the nucleus 
(per nucleon) and the hadron, respectively. The rapid- 
ity y and the transverse momentum q± of the lepton 
pair can be expressed in terms of the hadronic Man- 
delstam invariants S = (Pi + P 2 ) 2 , T = (Pi - q) 2 and 
U = (P 2 - q) 2 as: y = \\n{{Q 2 - U)/(Q 2 ~ T)) and 
ql = (Q 2 ~U)(Q 2 -T)/S~Q 2 . 

It is convenient to express the angular distribution of 
the lepton pair in terms of the helicity amplitudes W aiCr > 



of the virtual photon. A complete set of such amplitudes 
is R 



W T = Wi, u W L = W ,o, W AA = Wi,-! 
W A = -j= (Wi, + Wo.i) 



(4) 



which are referred to as the transverse, longitudinal, spin- 
flip, and double spin-flip structure functions. In terms of 
these helicity structure functions the cross section can be 
written as 

da a 2 



dQ 2 dq 2 ± dydn 64ir 3 SQ 2 



(wtl (1 + cos 2 6) + Wi Q - | cos 2 9 

+ Wa sin 29 cos + Waa sin 2 9 cos , (5) 

where Wtl = Wt + \ W\,. Integration over the angles 
reduces the cross section to 

da a 
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dQ 2 dq]_ dy 12tt 2 SQ- 
The ratio of (||) and (g) is often parametrized as 



(6) 



1 + A cos 2 9 + fi sin 29 cos <p + - sin 2 9 cos 2<j) . (7) 

We have verified that we reproduce the known results for 
the leading-twist matrix elements In particular, we 
find Waa = |Wd or X + 2u = 1, which is the Lam- Tung 
relation. 

For the calculation of the double-scattering contribu- 
tions we use the formalism of Luo, Qiu and Sterman [^|,^| . 
Here we present only our main results and refer the reader 
to a forthcoming publication for details of the calculation 
p"5| . We factorize the hadronic tensor for all contribu- 
tions in the usual way, and we keep only those matrix ele- 
ments, in which the two partons from the nucleus appear 
separately in combinations of color singlet operators. As- 
suming the absence of long-range color fluctuations in the 
nuclear wave function, only such correlators will show a 
nuclear enhancement. 

It is useful to distinguish two contributions to double 
scattering |@| ■ In the first case (double- hard process), 
both QCD interactions are hard and the beam parton 
can be considered as on-shell between the interactions. 
In the second case (soft-hard process) the beam parton 
picks up a soft nuclear parton before the hard scattering. 
The two contributions differ in the pole structure of the 
hard part of the scattering tensor. We neglect the inter- 
ference terms between soft and hard rescattering, as it is 
appropriate when q± is not too small. For q\ -C Q 2 the 
interference terms are important and eventually spoil the 
nuclear enhancement Q. 
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First we investigate the case of double-hard scatter- 
ing. The two-parton distribution functions depend on 
two Bjorken parameters, x a = Q 2 /{Q 2 — T) and Xh- The 
cross section is determined by three universal nuclear ma- 
trix elements that were introduced in || . We only denote 
the gluon-quark distribution function explicitly: 
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Here the superscripts ± denote the light-cone compo- 
nents of the four- vectors. The phase factors in (||) restrict 
the integration ranges in such a way that one integral over 
the longitudinal extension of the nucleus remains. This is 
the origin of nuclear enhancement. The quark- antiquark 
and two-quark distribution functions, X^ H (x a , Xh) and 
T^ i {x a , Xh), are given by similar expressions involving 
the correlator (q"f + qq^/ + q) ■ The two-gluon correlator 
does not contribute to double-hard scattering. 

The cross section is given by a convolution of a twist-2 
structure function / c /jj(0 of the hadron, the double hard 
matrix element and a perturbatively calculable hard part 
H ab+c (x a ,x h ,0- 
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(x a ,x h )H ab+c (x a ,x h ,0- (9) 



Here a, b and c denote the involved partons which have 
to be summed up. Note that Xh = Xh((.) depends on £. 

We hnd by explicit calculation that the Lam- Tung rela- 
tion holds for double-hard scattering. In the GJ frame we 
also hnd that the double-hard contributions do not gener- 
ate deviations from the simple (1 + cos 2 9) behaviour, i.e. 
fi = v = 0. The explanation for this surprising behavior 
is that the double-hard process resembles the classical 
picture of double scattering. For example, a quark from 
the hadron first scatters on a hard gluon from the nu- 
cleus, returns to its mass shell by radiating a gluon with 
large transverse momentum, and then annihilates on a 
hard antiquark from the nuclear quark sea. 

For soft-hard scattering one finds that the dominant 
contributions come from two universal two-parton dis- 
tribution functions, T^(x) and T^(x), which depend 
on a single Bjorken parameter x. They are obtained as 
integrals over the nuclear quark-gluon and two-gluon cor- 
relators, (F u +qj+qF + J and (F u+ F + U F X + F\) , respec- 
tively. Following Luo et al. H, we neglect the contribu- 
tions from two-quark correlators involving one soft quark, 
because these are far exceeded by the correlators involv- 
ing a soft gluon. Again, assuming the absence of long- 
range color fluctuations in the nucleus, an enhancement 
proportional to the nuclear radius is obtained. Explicit 
evaluation of the soft-hard cross section reveals that the 
Lam- Tung relation is violated by this contribution. 



Before we present numerical results, we need to spec- 
ify the nuclear twist-4 matrix elements. In the absence 
of experimental results, we follow p 16 and express the 
two-parton distribution functions T^^as products of the 
known one-parton distribution functions f(x): 



T^(x a ,x h ) = C(A)Af a/A (x a )f b/A (x h ) 



(10) 



Here f a / A is the parton distribution (neglecting shadow- 
ing and EMC-type effects) in a nucleus with mass number 
A. The normalization constant is assumed to grow like 
the nuclear radius, C(A) = CA 1 / 3 with C = 0.072 GeV 2 
[Q. For soft-hard scattering the twist-4 matrix elements 
only depend on the momentum of the hard parton. We 
assume that they are proportional to the distribution 
functions of the hard parton and include the effect of the 
presence of an additional soft gluon as a renormalization: 



(x) 



A^QSG^ 4/3 /a M (x) 



(11) 



where Alqsg = 0.1 GeV is taken from |17|. We use the 
twist-2 distribution functions from the CTEQ3M param- 
eterization for the proton [jl8| . 

We first discuss the application of our results to the 
data for 800 GeV protons on W from the experiment 
E866 at Fermilab |ll[] . Figure || shows the contribution 
from the single spin-flip amplitude Wa to the differential 
DY cross section (JsJ) in the CS frame. (Note that Wa van- 
ishes in p+p collisions for symmetry reasons.) The twist- 
4 contribution is much larger than the twist-2 result, and 
it is completely dominated by the double-hard scattering 
process. Note that the soft-hard contribution is negligi- 
ble. This feature can be tested by checking experimen- 
tally for a violation of the Lam- Tung relation which is 
exclusively given by soft-hard processes. Furthermore a 
characteristic property of our soft-hard and double-hard 
contributions is the specific difference between the CS 
and GJ frames 
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FIG. 2. The helicity amplitude w A = Q 2 W / a/(64tt 3 S'Q 2 ) 
for 800 GeV p + W at Q = 5 GeV and q± _ = 4 GeV. The plot 
shows the twist-2 result (solid line), the double-hard (dashed 
line) and soft-hard (dotted line) twist-4 contributions, in the 
CS frame. The soft-hard result is multiplied by a factor of 10. 
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would be very interesting to calculate the same effect in 
other approaches A measurement of the nuclear 

enhanced contribution should be possible at RHIC. 
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FIG. 3. The helicity amplitudes (a) totl and (b) wa 
in the CS frame for 250 GeV protons colliding with 
100 GeV/nucleon Au nuclei at Q = 5 GeV and q ± = 4 GeV. 
For details see Fig. |^. 

Next, we turn to p + Au collisions at the full RHIC en- 
ergy. As Fig. H shows, the twist-4 contribution is again 
larger than the leading- twist result. The twist-4 cross 
section peaks at negative rapidities (the direction of the 
proton beam), because the requirement of a high trans- 
verse momentum enhances the contribution, in which a 
quark from the proton annihilates with an antiquark from 
the Au nucleus. Again, the soft-hard process is negligible. 
The kinematics of the peak (Q — 5 GeV, q± = 4 GeV, 
y w — 2)fits well with the acceptance of the Phenix de- 
tector (g0[. The predicted count rate is not high (about 
30 events/year in a window of 1 GeV 4 ), but the char- 
acteristic forward-backward asymmetry should be easily 
detectable. 

In summary, we have calculated the nuclear enhance- 
ment of the twist-4 contribution to DY production in p + 
A collisions, with its full angular distribution. For RHIC 
energies the soft-hard contributions involving highly non- 
trival multi-field correlators are negligible, which is good 
news for any comprehensive QCD-description of high- 
energy heavy ion collisions. The dominant contribu- 
tion from double-hard scattering has an interpretation 
in terms of the classical double scattering picture. While 
double-hard scattering respects the Lam- Tung relation, 
soft-hard scattering does not. Double-hard scattering 
does not contribute in the GJ-frame, which allows for 
a sensitive experimental test of the LQS-formalism. It 
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